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The Pd(100)-(v / 5 X V5)R27 -0 surface oxide revisited 
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Combining high- resolution core-level spectroscopy (HRCLS), scanning tunneling microscopy 
(STM) and density-functional theory (DFT) calculations we reanalyze the Pd(100)-(v^> x VE)R27°- 
O surface oxide phase. We find that the prevalent structural model, a rumpled PdO(OOl) film 
suggested by previous low energy electron diffraction (LEED) work (M. Saidy et al, Surf. Sci. 
494, L799 (2001)), is incompatible with all three employed methods. Instead, we suggest the 
two-dimensional film to consist of a strained PdO(lOl) layer on top of Pd(100). LEED intensity 
calculations show that this model is compatible with the experimental data of Saidy et al. 
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I. INTRODUCTION 



Due to their importance for catalysis and corrosion, 
oxidation processes at transition metal (TM) surfaces 
have long received significant attention both in funda- 
mental and applied research. Roughly divided into pure 
on-surface (dissociative) adsorption, surface-oxide forma- 
tion, and oxide film growth, particularly the second step 
in the oxidation sequence is poorly understood on a mi- 
croscopic level. Not least, this is due to the fact that an 
atomic-scale investigation of TM surfaces on the verge of 
oxide formation still poses a significant challenge for the 
otherwise well-developed machinery of ultra-high vacuum 
(UHV) surface science: The rather high oxygen partial 
pressures and elevated temperatures required to initiate 
oxide nucleation, a low degree of order at the surface and 
highly complex, large unit-cell geometries name but a 
few of the problems encountered on the route towards a 
microscopic characterization of surface oxides. 

The vast structural parameter space connected with 
surface oxides is often prohibitive for an analysis based 
on one technique alone, in particular for the exhaus- 
tive searches required in diffraction-based structure de- 
termination techniques like low energy electron diffrac- 
tion (LEED). Fortunately, over the last 10 years it be- 
came clear that the joint effort of experimental works and 
theory (in particular when using first-principles methods) 
is synergetic and most valueable, if not crucial, and allows 
a convincing identification and characterization of com- 
plex, novel structures at surfacesi*2*Mi^. Such a multi- 
method approach is also adopted in the present work, 
where high- resolution core- level spectroscopy (HRCLS), 
scanning tunneling microscopy (STM), and density- 
functional theory (DFT) are employed addressing the 
Pd(100)-(V5 x y/5)R27°-0 surface oxide phase (coined 
\/5-phase in the following for brevity), which for T > 
400 K concludes the series of ordered phases observed on 
Pd(100) before three-dimensional cluster growth sets in&. 



A previous tensor LEED analysis suggested the >/E- 
phase to correspond essentially to a rumpled PdO(OOl) 
plane on top of Pd(100) 7 ^. In the following we will show 
that this assignment can not be reconciled with cither 
HRCLS, STM, or DFT, the methods employed in the 
present study. Instead, we propose a new model, consist- 
ing of a strained PdO(lOl) layer on Pd(100). This is in 
agreement with all experimental data and energetically 
more stable than the previous model. Moreover, we per- 
formed preliminary LEED calculations to show that this 
new arrangement can also be reconciled with the previous 
experimental LEED data of Saidy et a& Interestingly, 
the (101) orientation does not correspond to a preferred 
growth direction^ or a low-energy surface of crystalline 
PdO, suggesting that the film-substrate interaction may 
stabilize higher energy crystal faces. We argue that the 
atypical orientation and the resulting chemical proper- 
ties of such thin oxide layers may be of interest in future 
applications. 



II. EXPERIMENTAL AND COMPUTATIONAL 
DETAILS 

The STM measurements were performed in a UHV 
chamber with a base pressure below 1 x 10~ 10 mbar. The 
Pd(100) surface was cleaned by cycles of Ar + sputtering, 
annealing and oxygen treatments keeping the sample at 
900 K in an oxygen pressure of 2 x 10 -8 mbar followed by 
flashes to 1400 K. The cleanliness of the Pd(100) surface 
was checked by Auger Electron Spectroscopy (AES); no 
contaminants such as C and O could be observed within 
the detection limits. The -\/5-phase was thereafter formed 
by exposing the Pd(100) surface to an oxygen pressure 
of 5 x 10~ 6 mbar for 300 seconds at T = 600 K. 

The HRCLS measurements were conducted at the 
beam line I31tt2at MAXII in Lund, Sweden. The clean- 
ing procedure and preparation of the -\/5-phase was iden- 
tical to that described above. The HRCL spectra were 
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recorded at liquid nitrogen temperatures and at normal 
emission angle. The cleanliness of the Pd(100) surface 
was checked by monitoring the Pd 3d$/2, O Is, and C Is 
core-levels, as well as the valence band region; again, no 
contaminants could be detected. 

The DFT calculations were performed within the Full- 
Potential Linear Augmented Plane Wave (FP-LAPW) 
schemeiii^ii^ using the generalized gradient approxima- 
tion (GGA) 14 for the exchange-correlation functional. 
The y5 surface oxide was modeled in a supercell geome- 
try, employing a symmetric slab consisting of five layers 
Pd(100) in the middle plus the various PdO layers de- 
scribed below on both sides. A vacuum region of ~ 15 
A ensures the decoupling of the surfaces of consecutive 
slabs. All atomic positions within the PdO and the out- 
ermost substrate layer were fully relaxed. 

The FP-LAPW basis set parameters are as follows: 
i?P d T =1.8 bohr, i?° T =1.3 bohr, wave function expan- 
sion inside the muffin tins up to l^ ax = 12, potential 
expansion up to = 4, and a (4 x 4 x 1) Monkhorst- 
Pack grid with 8 k-points in the irreducible part of the 
Brillouin zone (28 k-points in the full zone). The energy 
cutoff for the plane wave representation in the interstitial 
region between the muffin tin spheres was -E™ f ax = 20 Ry 
for the wave functions and E™* = 169 Ry for the poten- 
tial. With this basis set very good agreement is obtained 
with the results of the structure determination for the 
p(2 x 2) phase on Pd(100) reported by a previous LEED 
study^, and we find the relative energetic stabilities of 
the various tested overlayer models converged to within 
±50meV per O atom. 

This does, however, not comprise the uncertainty in 
the absolute binding energies, Eb- DFT, even within 
the GGA, is known to poorly describe gas phase oxygen 
and gives in particular the binding energy for molecular 
O2 wrong by about 0.5 eV per O atom 16 . As this directly 
effects the obtained absolute binding energies for the var- 
ious surface oxide phases, a possible workaround would 
be to determine the total energy of molecular oxygen, 
Eq 1 , not via gas phase calculations, but viai£ 

l/2£#* « £^ Ojbnlk -£^ |bulk + AH / (300K s latm) 

i.e. employing an approximate equation for the PdO heat 
of formation, AHf, into which only the total energies of 
bulk PdO and Pd bulk, S*>° d t o bulk and i^*^, enter. 
Using the experimental A_ff/(300K, latm) one can thus 
arrive at £g ot without having to resort to atomic calcu- 
lations, though at the expense of discarding a completely 
first-principles type description. 

In the present work we will only compare the stability 
of various structural models all including the same num- 
ber of oxygen atoms. Then, the difference between the 
standard computation of binding energies, i.e. utilizing 
gas-phase computed Eq*, and the aforedescribed proce- 
dure amounts only to a constant shift in the calculated 
binding energies. Employing Ai7y Xp (300K, latm) = 
0.88 eV 18 , this shift amounts to 0.43 eV per O atom with 



a lower stability of the A7J/-derived binding energies and 
not including zero-point vibrations. This indicates the 
sizable uncertainty in the absolute binding energy val- 
ues and correspondingly dictates a cautious judgement 
on the endo- or exothcrmicity of a structure. In the fol- 
lowing, we will always indicate binding energies obtained 
with gas-phase computed Eq* (not including zero-point 
vibrations), which according to the above argument are 
likely to represent an upper limit to the real stability. Our 
main conclusion, the rebuttal of the prevalent v^-modcl, 
will, however, rather be based on relative energetic differ- 
ences, which are fortunately enough much better defined. 

The surface core-level shift (SCLS), Ascls, is defined 
as the difference in energy which is needed to remove a 
core electron either from a surface or a bulk atomic. In 
the initial-state approximation the SCLS arises simply 
from the variation of the computed orbital eigenenergies 
before the excitation of the core electron. In final-state 
calculations, on the other hand, the SCLS involves an 
additional component due to the screening contribution 
from the valence electrons in response to the created core 
hole, obtained approximately via the Slater- Janak transi- 
tion state approach of evaluating total energy differences 
using impurity type calculations as explained in detail in 
ref. |2(J. In case of the Pd 3d SCLSs the bulk level posi- 
tion can be employed as a well defined reference level to 
align theoretical and experimental spectra. For the Ols 
levels however, we note that it is not very practical to also 
use the Pd 3d bulk-level position to align the O Is data, 
as both types of orbital eigenenergies exhibit different 
convergence behavior with respect to the employed basis 
set. Also the Fermi-level position as another reference 
level present in both experimental and theoretical spec- 
tra is not very practical: Particularly for experiments on 
systems with high oxygen loads like surface oxides band 
bending can not be excluded. 

Fortunately, we only need to rely on the existence (or 
non-existence) of a split O Is spectrum in the present 
work. This difference in relative O Is level positions of 
various atoms within the same geometry is well defined 
and independent of the reference zero used. Just in or- 
der to present the theoretical and experimental data in 
the same plot (Fig. 2) we will therefore employ the sim- 
plest possible alignment approach given by equating the 
position of the lowest theoretical and experimental O Is 
core-level position. We stress that this crude procedure 
does not enter our physical argument and is solely used 
for graphical purposes. 



III. RESULTS AND DISCUSSION 

1. Shortcomings of the prevalent \ // 5-LEED model 

Fig. 1 shows the development of the experimental 
HRCL spectra from the O Is and Pd 3^/2 levels for the 
sequence of ordered structures that form on the Pd(100) 
surface with increasing oxygen coverage. After the two 
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FIG. 1: Experimental HRCL spectra of the O Is and the Pd 
3^5/2 levels for the sequence of ordered structures that form 
on Pd(100) with increasing oxygen coverage. The photon en- 
ergies were 650 eV and 400 eV respectively. 



known adsorption phases, p(2 x 2) at 9 = 0.25 ML and 
c(2 X 2) at 6 = 0.50 ML^i 21 i 22 i 2 ?i 24 , at first an interme- 
diate (5 x 5) surface oxide forms, before finally the v5 is 
obtained, on which we concentrate in this work. Focus- 
ing first on the Ols spectrum of this -\/5-phase, its most 
surprising feature is the existence of two sharp peaks in 
contrast to the single peak observed at all lower coverage 
structures. As is apparent from Fig. 1, a definite as- 
signment of these peaks is hampered by the energetically 
very close lying Pd 2>p levels. Still, we attribute these two 
peaks to emission from Ols levels as no similar features 
are observed at the low binding-energy side of the Pd 3d 
spectrum. Changing the incident photon energy to vary 
the escape depth of the photoelectrons leads us to con- 
clude that they originate from oxygen atoms close to the 
surface and at roughly the same depth. 

To arrive at a rough estimate of the relative coverages, 
we approximately remove the contribution from the Pd 
3p levels by subtracting the Pd 3p spectrum recorded 
from the clean surface. The corresponding HRCL spec- 
tra for the a/5- phase are shown in Fig. 2, now recorded 
at higher photon energies of 900 eV and 650 eV for the 
O Is and Pd 3d levels, respectively. At the expense of 
a decreased resolution, these somewhat high energies are 
employed to avoid diffraction effects and thus permit a 
rough estimate of the relative coverages giving rise to 
the two peaks. We note that despite the decreased res- 
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FIG. 2: Top panel: measured HRCL spectra of the O Is 
and the Pd 3d 5/2 levels from the Pd(100)-(v / 5 x y/5)R27°-0 
phase at higher photon energies (900 eV and 650 eV for O Is 
and Pd 3d, respectively) and with the Pd 3p contribution of 
the clean surface removed. Bottom panel: Calculated final- 
state shifts for the three structural models shown in Fig. 3. 
For Pd 3d the bulk-level is employed to align theoretical and 
experimental spectra. For O Is the lowest-energy theoretical 
peak is simply aligned to the lowest-energy experimental peak 
(see text). Note that only the PdO(lOl) layer on Pd(100) 
exhibits a split O Is spectrum with two significantly shifted 
components. See Fig. 3 for the nomenclature used to describe 
the atoms from which the various theoretical Pd core-level 
shifts originate. 



olution, the two peaks are still clearly distinguishable in 
the original (not shown) and subtracted (Fig. 2) data. In 
both cases we obtain a considerable binding energy shift 
of w 0.75 eV between both peaks and a ratio of about 
1:1 for the two components. On the basis of this analysis 
of the experimental O Is spectrum we would therefore 
anticipate at least two oxygen species at or close to the 
surface of the \/5-phase in close to equal amounts. 

Turning to the Pd 3d 5 / 2 spectra in Fig. [3 at least 
three oxygen-induced components are experimentally re- 
solved, at -0.32 eV, +0.38 eV and +1.30 eV, where a posi- 
tive SCLS indicates a higher binding energy with respect 
to the reference Pd bulk component. Using simple initial- 
state arguments, we expect an increased oxygen coordi- 
nation to yield a positive SCLS for the Pd 3c?5/2 level. 
Thus, the component shifted by +1.30 eV could be due 
to highly oxygen coordinated Pd atoms, agreeing with 
previous results from the intermediate oxide structure 
on Pd(lll) 5 , in which a very similarly shifted compo- 
nent was found and assigned to Pd atoms fourfold coor- 
dinated to oxygen. The component shifted by +0.38 eV 
could correspondingly be due to Pd atoms coordinated 
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PdO(001)/Pd(100) PdO(100)/Pd(100) PdO<I01)/Pd(lU0) 
E B = -0.12eV E„ = 0.90eV E B =l.35eV 
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FIG. 3: Top- and side-view of the three structural models 
of the Pd(100)-(V5 x VE)R27°-0 phase considered in the 
present work. All models assume an oxygen coverage of 0.8 
ML. a) PdO(OOl) layer on Pd(100)^, b) PdO(lOO) layer on 
Pd(100), and c) PdO(lOl) layer on Pd(100). The unit-cell 
is sketched in the top- views (solid line) , while the dashed lines 
indicate the direction of atomic rows seen in the STM images. 
The DFT binding energy of the three models clearly reveals 
the PdO(lOl) layer on Pd(100) as the most favorable model. 
2/ and 4/ denote two- and fourfold coordinated first-layer Pd 
atoms, and relate to the labels in Fig. 2 used to specify the 
atomic origin of the various computed core- level shifts. 



to less oxygen atoms (possibly two or three) , whereas the 
component shifted by -0.32 eV likely originates from Pd 
atoms at the interface between the Pd(100) substrate and 
the thin oxidic film. 

Trying to analyze the compatibility of the present 
HRCLS data with the prevalent structural model of the 
-\Z5-phase suggested on the basis of the tensor LEED 
analysis (henceforth abreviated with "LEED model" 
cf. Fig. 3a, we used the published atomic positions of 
the LEED model as input to our DFT computations. 
Obtaining SCLS that did not resemble the experimen- 
tal data at all, we initially proceeded by subjecting the 
LEED model to a complete structural relaxation. The 
resulting final-state SCLS after relaxation are shown in 
Fig. 2 and are still difficult to reconcile with the experi- 
mental data: The large splitting of the O Is spectrum is 
not reproduced and almost identical O Is positions are 
obtained for all O atoms in the structure. Recalling that 
the LEED model essentially corresponds to a PdO(OOl) 
overlayer on Pd(100), in which all oxygen atoms are in 
principle equivalent, cf. Fig. yp,, this result is not sur- 
prising. The agreement in case of the Pd 3d SCLSs is 
not much better, obtaining computed shifts that are split 
into two distinct groups in contrast to the four compo- 
nent structure seen experimentally. 

Geometrically, the structural relaxation in DFT al- 
most completely removes most of the strong rumpling 
introduced in the LEED study to fit the experimental 
/(-E)-curves, and the substrate/oxide interface smoothes 
out at a large interface distance indicating a very weak 
coupling. In the end, the absolute binding energy per 
O atom, Ebi of the LEED model was (with the caveat 
given in Section II) still found to be slightly endothermic 




FIG. 4: a) STM image showing a domain boundary between 
the Pd(100)-(\/5 x VE)R27°-0 and the p{2 x 2) phase, b) 
The same STM image but with a Pd(100) lattice superim- 
posed (grey circles). This shows directly that in the \/5-phase 
bright spots (assigned to Pd atoms, white circles) are shifted 
in neighboring rows by half a nearest neighbor distance, and 
that the dark spots (Pd atoms in hollow sites, transparent 
circles) coincide with hollow sites of the underlying Pd(100) 
substrate (Tunnel parameters: V=0.76 V, 1=0.57 nA). 

with respect to molecular oxygen {Eb = — 0.12 cV/ O 
atom), adding to our doubts of this prevalent structural 
model of the "\/5-phase. 

Finally, we also performed STM measurements of the 
\/5-phase, a corresponding image of which is shown in 
Fig. 3a. As may directly be seen from the image, neigh- 
boring bright rows arc shifted by half a nearest neighbor 
distance with respect to each other, in contrast to what 
would be expected from the PdO(OOl) geometry depicted 
in Fig. 3a. This finding (STM) together with the absence 
of a split in the calculated O Is spectrum (HRCLS) and 
the low energetic stability (DFT) led us to conclude that 
the prevalent LEED V5-model is incompatible with the 
three methods employed in the present study. 



2. Searching for a new model 

To identify an alternative geometry of the \/5-phase, 
we first further analyze our experimental STM data in or- 
der to reduce the vast phase space of possible structural 
models. The bottom right half of the STM image shown 
in Fig. 4a exhibits a domain of the coexisting p(2 x 2) 
on-surface adsorption phase, in which the oxygen over- 
layer simply occupies fourfold hollow sites of the underly- 
ing Pd(100) substrato 15 i 21 i 22 i 23 i 24 . Based on this known 
geometry and the frequent finding that oxygen would ap- 
pear as dark spots in the STM image, we construct the 
Pd(100) lattice and superimpose it on the experimental 
image, as shown in Fig. 4b. It then follows that the dark 
spots in the v^-phase are directly situated on top of the 
fourfold hollow sites of the Pd(100) lattice. 

Assuming that the bright protrusions in the STM im- 
age correspond to the geometric position of Pd atoms, we 
may further draw the Pd sublattice of the suspected sur- 
face oxide layer into the STM image as done in Fig. 4b. 
This way, a total of three Pd atoms per y/5 unit-cell are 
found, forming a rather open interlaced ring-like layer, 
the structure of which doesn't resemble a bulk-like PdO 
planar nearest-neighbor environment at all. The latter 
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would instead be obtained, if Pd atoms would also be 
present at the position of the large dark spots, yielding 
then a more compact layer with a total of four Pd atoms 
per y/E unit-cell. As the dark spots are directly situated 
on top of the hollow substrate sites, a straightforward ex- 
planation why these latter Pd atoms do not show up in 
the STM images would e.g. be a large corrugation within 
the surface oxide overlayer, in which all Pd atoms over 
Pd(100) hollow sites are strongly relaxed inwards. 

While we may thus tentatively determine the positions 
of the Pd atoms on the basis of the STM images, the lat- 
ter do not lead to any conclusions about the position 
and number of oxygen atoms in the y/E unit-cell. Con- 
cerning the O coverage, we can however resort to the 
HRCLS measurements. Calibrating the spectra, cf. Fig. 
1, with the p(2 x 2) (0.25 ML) and c(2 x 2) (0.50 ML) 
adsorbate structures known from the previous LEED 
worli 15 i 21 i 22 i 23 i 24 , a rough estimate of 6 ~ 0.8 ML is ob- 
tained, which would correspond to four O atoms per y/E 
unit-cell. 

Using these experimental observations, we proceed to 
set up structural models that are compatible with the 
data discussed so far. Assuming the v^-phase to be some 
form of surface oxide on Pd(100), PdO-like overlayers 
seem a most appealing choice for a model. A system- 
atic look at all possible low-index PdO planes turns up 
two orientations which exhibit Pd positions whose lateral 
arrangement would agree with that deduced from STM: 
PdO(lOO) and PdO(lOl), cf. Fig. 3b and 3c respectively. 
In contrast, the earlier LEED model consists essentially 
of a PdO(OOl) plane on Pd(100), which is not equivalent 
to PdO(lOO) due to the tetragonal unit-cell of PdO 25 . 
Hence, the LEED model features an orientation, which 
does not fit the STM data, as seen when comparing Fig. 
3a with Fig. 4. 



3. PdO(101)/Pd(100) as the new model 

Having filtered out PdO(lOO) and PdO(lOl) as two 
possible candidates for a new structural model, we note 
that the two differ only in the vertical position of the 
oxygen atoms. Each structure contains four oxygen and 
four Pd atoms per y/E unit-cell, which nicely fits the ex- 
perimental coverage estimate described in Section III. 2 . 
Yet, PdO(lOO) has all four O atoms above the Pd layer, 
while PdO(lOl) has two up and two down, as can be seen 
in Figs. 3b and 3c respectively. To discern between the 
two orientations, we subjected both to a full structure 
optimization in our DFT calculations. Interestingly, this 
yields a significantly increased stability for both over- 
layer models: in both cases, the binding energy is more 
than 1 eV per O atom higher than that of the previous 
LEED model, thus providing the final evidence for the 
incorrectness of the latter. More precisely, we find the 
PdO(101)/Pd(100) \/5-geometry to be the most stable 
of the three structural models depicted in Fig. 3 with a 
binding energy of Eb = +1.35 eV / O atom (compared to 




FIG. 5: Top- and side-view of the PdO(101)/Pd(100) model 
for the V5-phase on the basis of the DFT calculations. The 
rumpling of both the O and Pd atoms in the PdO overlayer 
and of the Pd atoms in the topmost substrate layer is given 
with respect to the center of mass of the respective layers. 
In the bottom right, also the average layer distances between 
these center of mass are indicated. 



E B = -0.12eV/O atom for the LEED model and E B 
= +0.90cV/O atom for the PdO(lOO) overlayer). 

Further checking on PdO(101)/Pd(100) as our new 
structural model, we also calculated a number of ge- 
ometries keeping the positions of the Pd atoms in the 
overlayer, but testing different lateral positions for the 
O atoms. In particular, this involved geometries, where 
some oxygen atoms were located in bridge sites between 
the overlayer Pd atoms to produce differently coordi- 
nated O atoms at the surface that could then possibly 
also generate a split O Is core-level spectrum. Yet, in 
all such combinatorial cases with O in hollow and bridge 
sites we obtained binding energies more than 0.2 eV lower 
than for our PdO(101)/Pd(100) model. As second test 
series, we tried different registries of the PdO(lOl) over- 
layer on Pd(100), i.e. we laterally shifted the PdO over- 
layer around on the substrate. Again, this always re- 
sulted in a lower stability with respect to the structure 
shown in Fig. 3c. Finally, the dark holes seen in the 
STM image suggested that one of the four overlayer Pd 
atoms could be absent, namely the one over the fourfold 
hollow site, Pdn in Fig. 5. With bulk Pd as reservoir for 
the removed Pd atom, we also find such a model to be 
less stable. In conclusion, the DFT calculations there- 
fore strongly favor the PdO(101)/Pd(100) structure for 
the -s/5-phase. 

To see whether this new model is also compatible with 
the presented experimental HRCLS data, we show its 
computed final-state SCLS in Fig. 2 and list all its Pd 3d 
shifts in Table I. Comparing with the shifts obtained for 
the LEED model and for PdO(100)/Pd(100), cf. Figs. 
3a and 3b, PdO(lOl) is the only model that exhibits 
an appreciably split Ols core-level spectrum, due to the 
presence of both on- and sub-surface O in the geome- 
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Initial Screening Final 


Experiment 


I4/ , Pd Br /H 


+0.79 +0.10 +0.89 
+0.86 +0.05 +0.91 


+1.30 


hf, PdBr 

Pd H 


+0.20 +0.09 +0.29 
+0.18 +0.21 +0.39 


+0.38 


II, Pdl 


-0.31 +0.18 -0.13 
-0.28 +0.19 -0.09 
-0.26 +0.15 -0.11 
+0.01 -0.05 -0.04 
+0.16 -0.11 +0.05 


-0.32 



TABLE I: Calculated and measured Pd 3d surface core-level 
shifts for the PdO(101)/Pd(100) model in eV. The computed 
values are separated into initial-state and screening contribu- 
tion, yielding the total final-state shift that can be compared 
to experiment. See Figs. 2 and 5 for the notation to describe 
the various first (I) and second (II) layer atoms. 



try, cf. Fig. 2. The obtained Pd 3d shifts of +0.9 eV 
and +0.4 eV due to fourfold and twofold oxygen coordi- 
nated Pd atoms in the PdO(lOl) overlayer, Pde r /H an d 
Pder/PdH hi Fig. 5, compare reasonably with the two 
experimental peaks that had already been assigned to 
differently coordinated Pd atoms on the basis of initial- 
state arguments. The remaining experimentally resolved 
peak with a small negative shift had similarly been at- 
tributed to the top Pd substrate atoms at the interface, 
which in the calculations exhibit almost vanishing SCLSs 
(Pdl in Fig. 5). 

Of course, the very structure of the PdO(lOl) over- 
layer with an equal amount of on- and sub-surface oxy- 
gen atoms renders the measured splitted O Is spectrum 
immediately plausible. In fact, the significant rumpling 
together with the different sub-surface O coordination to 
the underlying substrate, cf. Fig. 5, yield even slightly 
different shifts for the two atoms of each oxygen species 
present in the -\/5 unit-cell, cf. Fig. 2. Averaging the 
contributions within each group, we obtain a computed 
initial (final) state shift of 0.55 eV (0.49 eV) between the 
O Is peaks due to on- and sub-surface oxygen atoms, 
in reasonable agreement with the measured value of 0.75 
eV. 



of experimental LEED data published by Saidy et al& 
(scanned and digitized from their Fig. 2). From this we 
may judge whether or not our geometry can yield LEED 
1(E) spectra on par or superior to those of Saidy et al. 
Our quantitative LEED calculations utilized the TensEr- 
LEED program package^, employing 10 fully relativistic 
phase shifts and a first-principles, energy-dependent real 
part of the inner potential, both generated for the sur- 
face geometry of Fig. fusing Rundgren's phase shift pro- 
gram packaged. Where the Tensor LEED method 28 i 29 
was employed, care was taken to ensure the full-dynamic 
reproducibility of the results in the final step of the calcu- 
lation. All non-structural parameters of the calculation 
were kept fixed at the values chosen by Saidy et al. 

As a first step, we simply used the exact optimized ge- 
ometry of the DFT-GGA calculations as input to the 
full-dynamic part of the TensErLEED code. Already 
this produced 1(E) curves in reassuring visual agreement 
with the scanned experimental spectra - i.e., all major 
spectral features could be reproduced. Still, shifts be- 
tween individual peaks and overall shape difference only 
allowed for an average Pendry R- factor 30 Rp = 0.51 be- 
tween calculated and scanned 1(E). Hence, in a sec- 
ond step, we used the Tensor LEED method to relax 
all vertical positions in the PdO(lOl) layer, as well as 
the topmost two Pd(100) substrate layers below. The 
result of this is a clear drop of the best-fit Pendry R- 
factor to Rp = 0.28, shared by both integer (i?i„t = 
0.29) and fractional (i?f rac = 0.28) beams on average. 
The improvement is mainly due to overall slightly ex- 
panded distances between the individual layers compared 
to the DFT-GGA result. Moreover, a significant buck- 
ling is found in the second substrate layer, which was 
not relaxed in DFT. Of course, some differences of this 
kind must be expected already because no lateral or non- 
structural degrees of freedom were adjusted in the LEED 
fit. While the latter parameters may well account also 
for the remaining discrepancies between calculated and 
experimental 1(E) curves, the main goal of our LEED 
calculations has clearly been achieved: Already a very 
limited structural refinement of PdO(101)/Pd(100) pro- 
duces experiment-theory agreement at a level which is 
even slightly improved compared to that presented by 
Saidy et al. (Rp — 0.306)& in their analysis. The consis- 
tency of our model with all available experimental data 
is thus established. 



4- Compatibility with existing LEED data 

So far, we have shown that our new model is superior 
to the model of Saidy et ali&> with respect to all three 
techniques employed in the present work (HRCLS, STM, 
DFT). On the other hand, the model of Saidy et al. has 
strong backing from quantitative LEED. Hence, a final 
verification of our model would be to establish its viabil- 
ity also by this method. To accomplish this goal, we per- 
formed a rather restricted set of LEED 1(E) calculations 
for our structural model, comparing them to the very set 



5. Strained PdO(101)/Pd(100) 

The new structural model for the v^-phase is essen- 
tially a strained and rumpled PdO(lOl) film on top of 
Pd(100). The PdO(lOl) in-plane lattice constant is al- 
most equal to that of a \/5 unit-cell on Pd(100), with 
the unit surface area of the commensurable film found 
here smaller by only 1.4% than for unstrained PdO(lOl). 
On the other hand, we compute a rather strong coupling 
of 100 meV/A 2 of the laterally compressed PdO over- 
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layer to the underlying Pd(100) substrate, rationalizing 
the formation of a commensurable surface oxide struc- 
ture. This strong coupling also helps to stabilize the 
particular PdO(lOl) orientation, which is experimentally 
not found to be a preferred growth direction of PdO 
crystallites^. Our calculations show that the stoichiomet- 
ric termination of bulk PdO(101)-(l x 1) suggested in the 
\/5-film, i.e. the one terminated by O atoms as shown 
in Fig. 5, is in fact considerably more stable than the 
two other ways of truncating PdO in (101) direction, 57 
meV/A 2 compared to 134 meV/A 2 (also O terminated) 
and 128 meV/A 2 (Pd terminated)^. Interestingly, the 
bulk PdO(lOO) orientation shown in Fig. 3b exhibits even 
a significantly lower surface energy (33 meV/A 2 ), while 
this orientation is in the commensurable thin film geom- 
etry discussed here energetically not as favorable as the 
PdO(lOl) v^-model. Evidently, the presence of oxygen 
at the oxide/metal interface yields a stronger coupling 
to the underlying substrate and is ultimately responsible 
for the higher stability of the PdO(101)/Pd(100) surface 
oxide geometry. 

This example of the stabilisation of a higher energy 
crystal face in thin oxide films due to strong interfacial 
coupling to the substrate adds another interesting aspect 
to the new physics found recently in studies concerning 
oxide formation at TM surfaces. Among other findings, 
the formation of incommensurable domains of low en- 
ergy oxide faces has been reported for ruthenium single 
crystals^ 2 ^, delineating the opposite case to the results 
reported here, i.e. when the oxide orientation is more im- 
portant than a good coupling to the underlying substrate. 
Apparently, the lower thermal stability of palladium ox- 
ides compared to Ru02 increases the importance of the 
oxide/metal interface. This is further supported by the 
surface oxide structure just found on Pd(lll), which does 
not resemble any PdO bulk orientation at al£. 

Experimentally, oxide thicknesses below about 20 A 
have been found in all of these cases, indicating either a 
slow growth kinetics once the thin films have formed or 
a thermodynamic hindrance to form thick bulk oxides. 
This could be of interest in oxidation catalysis, where 
such oxide patches forming on TM surfaces in the reac- 
tive environment are now discussed as the actually ac- 
tive material22iMiSiSi. If a continued growth of these 
oxide films is not possible, so that their structure al- 
ways remains significantly affected by the interfacial cou- 
pling, they may exhibit catalytic behavior which is non- 
scalable from corresponding bulk oxide crystallites - in 
other words, truly nano-catalytic properties. Even when 
the oxide growth is not limited, the structure of the ini- 
tially formed oxide film will be crucial, setting the stage 
for the ensuing oxidation process. For thicker films, in- 



terfacial coupling will be progressively less influential, so 
that an initially stabilized higher energy oxide orientation 
as found in the present work should eventually become 
liable to faceting. The corresponding three-dimensional 
cluster growth has indeed been observed for the contin- 
ued oxidation of both Pd(100) and Pd(lll)£i2£. 

IV. SUMMARY 

In conclusion, we have shown that the prevalent struc- 
tural model for the Pd(100)-(%/5 x \/5)i?27 o -0 surface 
oxide can not be reconciled with neither the experimen- 
tal nor the theoretical methods employed in the present 
study: Its surface symmetry does not fit to the one ob- 
served by STM, and the calculated HRCLS for this struc- 
ture do not show the appreciable splitting of the O Is 
spectrum observed experimentally. In addition DFT cal- 
culations give only a very low energetic stability and a 
relaxed geometry that does no longer exhibit the signif- 
icant rumpling originally introduced to match the mea- 
sured LEED 1(E) curves. 

Based on the present experimental data we reana- 
lyze the -s/5-phase and suggest an alternative structural 
model: a strained PdO(lOl) layer on Pd(100). This ar- 
rangement is energetically much more stable in our DFT 
calculations. Its computed final-state SCLSs agree well 
with all HRCLS measurements, linking the large split- 
ting of the O Is spectrum to the presence of oxygen both 
at the surface and at the oxide/metal interface. Already 
a very restricted set of LEED intensity calculations es- 
tablishes the compatibility of this structure also with the 
previously published LEED intensity data. 

The PdO(lOl) orientation, which is experimentally not 
found to be a preferred PdO growth direction, is stabi- 
lized by the strong coupling to the underlying substrate 
in the present thin film limit. In comparison to the en- 
suing three-dimensional cluster growth during continued 
oxidation, the -s/5-phase is therefore likely to display dif- 
ferent physico-chemical properties, which might be of 
interest or relevance to high-pressure applications like 
catalysis. 
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